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Based on assumptions which have led to the heat agreement between 
theory and test data on inelastic buckling of flat plates, a general set 
of equilibrium differential equations for the plastic buckling of cylin- 
ders has been derived. These equations have been used to obtain solu- 
tions for the compressive and torsional buckling of long cylinders in 
the yield region. 

Test data are presented which indicate satisfactory agreement with ■ 
the theoretical plasticity-reduction factors in most cases. Where a 
difference in results exists, test data are in substantially better agree- 
ment with the results obtained by use of the maximum-shear law rather 
than the octahedral- shear law to transform axial Btress-strain data to 
shear stress-strain data. 


INTROHJCTION 

Inelastic Compressive Buckling of Flat Plates 


The state of knowledge up to 1936 concerning inelastic buckling of 
plates and shells has been summarized by Timoshenko in reference 1 . The 
main efforts were concerned with attempts to modify the various bending- 
moment terms of the equilibrium differential equations by the use of 
suitable plasticity coefficients determined from experimental data on 
columns. Although such semiempirical efforts met with a reasonable degree 
of success, the theoretical determination of plasticity-reduction factors 
for flat plates has been achieved within recent years as the result of 
the development of inelastic-buckling theory. Because such developments 
are recent and form the background for the inelastic-buckling theory for 
nl-iftl 1 b developed herein, the following discussion concerning the assump- 
tions and results of the various theories is presented in some detail. 

Different investigators have used differing assumptions in the 
development of their theories. The major assumptions underlying each 
of these theories are given in the following table. 
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Investigator 

Stress- strain law 

Plasticity law 

Buckling model 

Bijlaard 
(ref. 2 ) 

Incremental and 
deformation types, 
v instantaneous 

Octahedral 

shear 

No strain 
reversal 

Ilyushin 
(ref. 3 ) 

Deformation type, 
" v = 0.5 

Octahedral 

shear 

Strain 

reversal 

Handelman and Prager 
(ref. 4) 

Incremental type, 
v instantaneous 

Octahedral 

shear 

Strain 

reversal 

Stowell 

(refs. 5 and 6 ) 

Deformation type, 
v = 0.5 

Octahedral 

shear 

No strain 
reversal 


Historically, Bijlaard (ref. 2) appears to have been the first to 
arrive at satisfactory theoretical solutions for inelastic -buckling 
theories. His work is the most comprehensive of an those considered 
in that he considers "both, incremental and deformation theories and con- 
cludes that the deformation type is correct since it leads to lower 
inelastic "buckling loads than those obtained fpam incremental theories. 
His work was first published in 1937* That paper and later publications 
include solutions to many important inelastic-buckling problems. How- 
ever, this work appears to have remained unknown to most of the later 
investigators . 

Ilyushin (ref. 3 ) briefly referred to Bijlaard' s work and then pro- 
ceeded to derive the basic differential equation for inelastic buckling 
of flat plates according to the strain-reversal model. The derivation 
of this equation is rather elegant and was used by Stowell (ref. 5), who, 
however, used the no-strain-reversal model. The differential equation 
obtained by Bijlaard reduces to that derived by Stowell by Bettin g 
v = l/2 in the former. Handelman and Prager (ref. 4), during this time, 
obtained solutions to several inelastic-buckling problems by use of 
incremental theory. Test data on compressed flanges and plates indi cate 
that the results of incremental theories are definitely unconservative 
regardless of the buckling model, whereas deformation-type theories are 
in relatively good agreement. 

The problem of plastic buckling has also been the subject of much 
experimental research. The use of the secant-modulus-reduction factor 
was first proposed for plates under c omp ressive loads by Gerard (ref . J) 
on the basis of tests on Z- and channel sections. Later, Stowell (ref. 5 ) 
proved theoretically that use of the secant modulus is correct for hinged 
flanges and that for elastically restrained flanges and plates the 
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plasticity-reduction factor Includes a function of the tangent modulus 
in addition to the secant modulus. For long columns, the factor depends 
only upon the tangent modulus. The plasticity-reduction factors proposed 
hy St owe 11 for simply supported flanges and plates have received excellent 
experimental confirmation (refs. 8 to 10) and it has "been well known for 
some 50 years that the tangent modulus is in good agreement with test data 
for col umns . Thus, the theoretical plastic-hunki ing factors for plates 
under compressive loads appear to he well substantiated hy rather precise 
experimental data. 


Inelastic Shear Buck! ing of Flat Plates 

In contrast with plastic compressive buckling, shear buckling of 
plates appears to he on less substantial ground. Ab the result of a 
series of tests on long 2024-0 aluminum-alloy plates under shear, Gerard 
(ref. 11) proposed use of the shear secant modulus as the plasticity- 
reduction factor for this case. The shear secant modulus is determined 
from a shear stress-strain curve, which, according to reference 11, is 
to he derived from an axial stress-strain curve on the basis of the 
maximum- shear plasticity law. St owe 11 (ref. 6) derived a theoretical 
plasticity-reduction factor for shear which has virtually the same numeri- 
cal value for all conditions of elastic restraint between simple support 
and clamped. 

In comparing the test data of reference 11 with the theoretical 
reduction factor, Stowe 11 used a shear stress -strain curve derived hy 
the octahedral-shear plasticity law. The shear plastic-buckling test 
data were found to lie consistently below the theoretical factor. Fur- 
thermore, St owe 11 attempted to explain the agreement between the shear 
secant-modulus method proposed in reference 11 and the test data therein 
on the basis that the stress -strain curve for 2024-0 aluminum alloy can 
be well approximated by a power law. 

Recently, in a series of tests on long, square, 2014-T6 aluminum- 
alloy tubes in torsion, Peters (ref. 10) presented a new set of test 
data on plastic shear buckling. Although the stress-strain curve of 
this material c anno t be adequately approximated by a power law, excellent 
agreement was found to exist between the new test data and the shear 
secant -modulus method proposed’ in reference 11. The theoretical factors 
of Stowell (ref. 6) and Bijlaard (ref. 2) were found to be consistently 
higher than the test data by an order of approximately 15 percent in the 
buckling stress. 

In summarization, then, the assumptions which lead to the best agree- 
ment between theory and test data on inelastic buckling of aluminum- alloy 
flat plates under compression loading are deformation-type stress-strain 
laws, stress and strain intensities defined by the octahedral-shear law. 
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and the no-strain-reversal model of inelastic "buckling. Although there 
may he theoretical objections to deformation theories as a class and the 
use of a no-strain-reversal model in conjunction with classical stability 
concepts, test data do suggest the use of results obtained from a theory 
based on these assumptions . 

For the inelastic buckling of flat plates under shear loading, 
plastic-buckling theory and the test data are not in good agreement . 

The principal difficulty appears to lie in the use of the octahedral- 
shear law to transform stress-strain data under axial loading to shear 
stress-strain data. This situation is discussed further herein in con- 
nection with results obtained for torsional buckling of cylinders. 


Inelastic Buckling of Cylinders 

Timoshenko (ref. 1) has presented some attempts to describe the 
inelastic buckling of a cylinder under axial ccmpreBsive forces subjected 
to axisymmetric buckling. These results are based on the intuitive use 
of the reduced modulus in place of the elastic modulus where the latter 
appears in the elastic-buckling-stress equation. 

Bijlaard (ref. 12), in an extension of his theory for inelastic 
buckling of flat plates, has considered the inelastic buckling of a 
cylinder subject to compression. Both the axisymmetric and the circum- 
ferential modes of buckling were considered in this analysis. The 
results are discussed in subsequent sections of this report. 

In this paper, a general set of equilibrium differential equations 
for the plastic buckling of cylinders is derived. This set of equations 
is perfectly general and applies to any loading system leading to buckling. 
In particular, solutions are obtained for compressive and torsional buck- 
ling of long cylinders in the yield region. 

The plasticity terms appearing in the equilibrium equations depend 
upon the choice of the buckling model. For the no-strain-reversal model, 
which is used in this analysis, the fact that the axial load must increase 
slightly during buckling in order that no unloading . should occur presents 
a mathematical difficulty when using classical stability concepts in which 
the loading remains constant during buckling. This difficulty is dis- 
cussed and an attempt to remove it is presented. 

This investigation was conducted at New York University under the 
sponsorship and with the financial assistance of the National Advisory 
C ommit tee for Aeronautics. 



NACA TO 3726 


5 


SYMBOLS 


plasticity coefficients defined by equations (A13) 


axial rigidity, E s t j ^1 - v 2 ) 


plate width 


bending rigidity, E s t^l2^1 - v 2 ) 


diameter 


modulus of elasticity 


secant modulus 


tangent modulus 


strain intensity defined by equation (A2) 


f x >Wz force 


shear elastic modulus 


shear secant modulus 


shear buckling coefficient 

length of cylinder 

bending moment per unit width 

number of longitudinal half wave lengths 

loading per unit width 

number of circumferential wave lengths 

external pressure 

radius of cylinder 


thickness 
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u,v,w displacements 

x, 8 , z coordinate s 

a = (3/<Xi 2 ) [l - (Et/E s )J 
P = t 2 /l2R 2 

7 shear strain 

e axial strain 

t] plasticity-reduction factor 

Tv = l/m 

v Poisson's ratio 

v e elastic value of Poisson's ratio, equal to 0.3 

a axial stress 

0j_ stress intensity defined "by equation (Al) 

t shear stress 

<P = Tcrl 1 " v2 )/ e s 
X curvature 

operator, j^d^/dx 2 ) + (^ 2 /r 2 S9 2 JJ 

V s - {Jf 

( )' variations which arise during buckling, such as M' and N’ 

Subscripts: 
c compression 

cr critical 


e 


elastic 
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f falling 

s shear 

x,y,xy coordinate orientation for M, N, a, and € 

1,2,3 variations which arise during "buckling in e and X 


EQUILIBRIUM EQUATIONS 


The plan followed in this report is to present the theoretical deri- 
vations in appendixes. The theoretical results, comparisons with test 
data, and a discussion of the significance of these results appear in 
the main body of the report. 

In appendix A, the assumptions of the plasticity theory used are 
discussed and the stress and strain intensities are defined according 
to the octahedral-shear law. Considerations involved in the buckling 
model are then considered and the incr eme ntal forces and moments which 
arise during buckling are presented based on the no-strain-re versa! model. 

In appendix B, Donnell 1 s (ref. 13 ) simplified strain-displacement 
and equilibrium equations derived originally for cylinder elastic-buckling 
problems are combined with the incremental force and moment relations of 
appendix A. In this manner, a complete set of equilibrium differential 
equations is obtained for use in the solution of cylinder plastic-buckling 
problems. Included in appendix B is an attempt to remove the difficulty 
of using equilibrium equations based on classical stability concepts for 
inelastic-buckling problems in which the no-Btrain-re versa! model requires 
that the load must increase slightly during buckling. 


AXIAL COMPRESSIVE BUCKLING OF A LONG CYUNDER 
Solution of Problem 


In appendix C, the axi symmetric buckling of a long cylinder under 
axi a l compression is considered. The critical stress obtained by use of 
the equilibrium equations derived in appendix B has the following form: 

®cr - IK 1 - ^)]" l/aE s( E fc/ E S ) l/2t / K 


( 1 ) 
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As in the inelastic buckling of flat plates, an effects of exceeding 
the proportional limit are incorporated in a plasticity-reduction factor 
defined as follows: 


% “ a cr/( ff cr) ( 


( 2 ) 


The solution for the elastic case is obtained by substituting 
E s = E, Et/E s = 1 , and v = ve in equation (l): 

( ff cr) e = [4 - v e 2 )J ^Et/R = 0 . 6 Et/R 

By use of equation ( 2 ), 


( 3 ) 


In general. 



(*) 


ffcr - O.firj^CfcAO 


( 5 ) 


Bijlaard (ref. 12) has previously obtained results equivalent to 
equations (l) and (4). His results are more exact, in fact, since the 
variation of Poisson's ratio in the inelastic range is included directly 
in the analysis. In the interests of simplicity, the. present analysis 
utilizes the artificial device of taking v = l/2 in both the elastic 
and plastic regions and then employs an approximate correction (see 
eq. (C15)) which, yields exact solutions for the elastic and plastic 
ranges as limits . This method follows a suggestion of Stowell (ref. 5 ). 

The purpose here in obtaining this solution was to present a unified 
approach to the campresBive and torsional buckling of cylinders from the 
set of equilibrium equations derived in appendix B. Furthermore, equa- 
tions (l) and (li-) serve as a basis for interpreting test data presented 
herein on the inelastic buckling of c ylinde rs. 

In addition to the axisymmetric case considered here, Bi^laard 
(ref. 12) has also considered the circumferential buckling mode of a 
cylinder under axial compression. As in the elastic case, the critical 
stresses obtained for both buckling modes are essentially the same. 
Bijlaard has pointed out that for mild steel the circumferential mode 
may lead to a slightly higher buckling load than the axisymmetric mode. 
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It would be a relatively simple matter to solve the circumferential 
buckling case by use of the equilibrium equations of appendix B. However, 
in view of. Bijlaard's results and the fact that the axi symmetric mode is 
often observed in tests on round cylinders which, buckle plastically, the 
solution for this buckling mode was not pursued. 


Test Data 

Osgood (ref. l4) and Moore and Holt (ref. 15) have presented test 
data on the falling strength of drawn circular tubes under axial com- 
pression. Osgood tested 2017-T^ aluminum-alloy tubes for which compres- 
sive stress-strain curves were given and also chrome -molybdenum tubes 
for which, unfortunately, neither the compressive stress-strain curves 
nor compressive yield stresses were obtained. It is interesting to note 
that photographs of the test specimens indicate the appearance of axisym- 
metrie buckling in some cases. Moore and Holt tested 6o6l-T6 aluminum- 
all oy tubes for which the compressive yield stresses were given although 
not the compressive stress-strain curves. A typical compressive stress- 
strain curve for this material, with a corresponding yield stress, was 
taken from reference 1 6 for correlation purposes. 

To reduce the experimental data for comparison with theory, the 
experimental failing strength was divided by the critical elastic stress 
(eq. (3)) to determine the experimental plasticity-reduction factors 
given in table 1. The theoretical value of tj c was determined in each 

case by use of the pertinent stress-strain data according to equation (4). 
For Poisson's ratio, the following values were used: 

v e = 0.3 (6) 

v = 0.5 - (E b /e)( 0.5 - v e ) 

The relation for the variation in Poisson ' b ratio in the yield region 
has been shown to apply to isotropic, plastically incompressible solids 
by Gerard and Wildhora (ref. 17). 

The theoretical and experimental values of q c as a function of 

the inelastic-compressive -buckling stress are plotted in figure 1. For 
the limited range of test data on 2017-T4 and 6o6l-T6 aluminum-alloy 
tubes it can be observed that good agreement is obtained. This is remark- 
able in view of the fact that while it tends to confirm the relation for 
the plasticity-reduction factor (eq. (if)) it also confirms the classical 
small-deflection stability theory for compressed cylinders. 
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TORSIONAL BUCKLING OF A LONG CYLINDER 
Cylinders of Extreme Length. 


In appendix D, the torsional buckling of a long cylinder is con- 
sidered. If the cylinder is of sufficient length, the boundary condi- 
tions at the ends have negligible influence and two-lobe buckling occurs. 
The critical stress for this case obtained by use of the equilibrium 
equations of appendix B has the following form: 

x _ 3/4 3/2 

Tcr = 0.272(l - v 2 ) E B (t/R) (7) 

In the elastic case, equation (7) becomes 

( T cr) e = 0.272(l - v e 2 )‘ 5/4 E(t/R) 5 ^ (8) 

Donnell (ref. 13) has shown that equation (8) applies for 


am 


where a = k 2 for simply supported ends and a = 60 for clamped ends. * 
By use of the equivalent of equation (2) for the torsional buckling case 


A-v e af A , 


In general, therefore. 


Tcr = 0.272(l - v e 2 )~^\^t/R )^ 2 


The plasticity-reduction factor depends primarily upon the secant 
modulus, which has been found to be the case whenever buckling occurs as 
a twisting action. This has been previously observed for compressive 
buckling of hinged flanges where buckling occurs as a twisting action 
and the plasticity-reduction factor depends primarily upon the secant 
modulus . 
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In appendix A, the secant modulus Is defined as 

E s = a i/ e i (12) 

If the shear stress- strain curve is known (this can he relatively simply 
obtained by torsion tests on tubes although no such direct ^experimental 
procedure exists for obtaining these data for flat plates), then the 
respective stress and strain intensities according to equations (Al) 
and (A2) are as follows: 

di = (3 ) 1//2 t ' (130 

®i = (3)“ 1 / 2 7 (l4) 

By use of equation (12) 

E s = 3t /y (15) 

Since -x/y = G fl , from equation (15) 

Gb = Es/3 

Therefore, equation (10) can be interpreted as 

s-f&U 

in cases in which the shear stress-strain curve is available. 

If, on the other hand, shear stress-strain data may not be avail- 
able, then it is possible to construct a shear stress-strain curve from 
simple tension and compression stress-strain data by use of either the 
octahedral- or the maximum-shear law. As indicated in reference 11 for 
materials which are anisotropic as a result of straightening (Bauschinger 
effect), it is probably best to utilize an axial stress-strain curve 
which is an average of the tension and compression curves each at 45° 
with the direction of the applied shear. If the axial stress in the 
simple axial test is 0 ^, then the stress and strain intensities according 
to the octahedral-shear and maximum-shear laws are 


(16) 


( 17 ) 


12 


MCA 03)1 3726 


°i " °x 
e i = e x 


(18) 


Thus, in "both cases. 


E s - °x/ e x (!9) 

However, the shear stresses corresponding to a x are different. 
For octahedral shear 


T = (3)'^% (20) 

and for maximum shear 

t = a x /2 (21) 

Thus, the values of E B and therefore ti s for a given value of a x 
correspond to a lower value of t for the maximum-shear law as compared 
with that of the octahedral-shear law. 

Cylinders of Moderate Length 

A solution for the inelastic buckling of cylinders of moderate 
length, in which case the boundary conditions have a decided influence 
upon the buckling stress, has not been obtained herein. However, the 
elastic solution is known and has been given in the following form by 
Batdorf (ref. 18): 


T cr = 0 . 7V7E (t /r (r/ 1 ) 1 ^ 2 


( 22 ) 


for 50t/R < (if R) 2 < lGR/t. 

For short cylinders defined approximately by 2^/nt < 1, the flat- 
plate solution 


T er - 




12 


(Wf 


(23) 
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applies where 


tl fl = Gs/G (24) 

The plasticity-reduction factor for this case has been proposed by Gerard 
(ref. 11 ) and is based on the maximum- shear law to transform the axial 
stress-strain data to shear data- 

Since equation (24) for short cylinders depends primarily upon the 
secant modulus as does equation ( 10 ) or ( 17 ) for long cylinders, it 
appears that equation ( 10 ) or (if) may be used as a reasonable approxi- 
mation of t] b for cylinders of moderate length. 


.Test Data 

Stang, Ramberg, and Back (ref. 19), Moore and Paul (ref. 20), and 
Moore and Bolt (ref. 15) have presented test data on the torsional failing 
strength of long and moderate -length drawn circular tubes. In most cases 
failure occurred as a result of inelastic buckling in the two-lobe mode. 

Stang, Bamberg, and Back tested 2017-T4 and chrome -molybdenum tubes 
for which representative shear s tre s s - strain curves were presented. In 
this case, therefore, it was possible to correlate theory and experiment 
on the basis of equation (17) . Of the large mass of test data given in 
reference 19 , a relatively small amount was useful for correlation pur- 
poses. These data are given in table 2 and were selected on the basis 
that the tubes were long in the sense of equation ( 9 ) (clamped ends). 
Furthermore, the yield stress of the 2017-T4 tubes was approximately 
25 ksi and the failing stress was less than 2 6 ksi to correspond with 
the given shear stress-strain data. For the chrome -molybdenum tubes, 
the corresponding values were 49 and 58 ksi, respectively. Many of the 
other test data were beyond the range of the given stress-strain data. 

In computing the experimental values of t} s , the torsional failing 

stress was divided by the critical elastic shear stress of equation ( 8 ). 
Since the shear stress-strain data were given, the theoretical values 
of tj b were calculated by use of equations ( 17 ) and ( 6 ) for Poisson’s 

ratio with G B jG replacing E B ^E in the latter. 

The theoretical and experimental values of t] s are shown in fig- 
ure 2. It can be observed that good agreement is obtained for the 
2017-T4 data, whereas the agreement is not so good for the chrome - 
molybdenum-tube data. 
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Moore and Paul (ref. 20) tested 605I-T6 seamless tubes of moderate 
length. Tension stress-strain data are also given. The torsional failing 
stresses are listed in table 3 together with the experimental value of T) s 

calculated by use of equation (22) for tubes of moderate length. The 
theoretical values of tj s were determined by use of equations (10) 

and (6), using the tension stress-strain data. Shown in figure $ la a 
comparison of the t} s values based on the use of the octahedral- and 

maximum-shear laws to transform the axial stress-strain data. Better 
correlation is obtained with the latter, although the test data are below 
the theoretical values of q s . Thi s may possibly reflect the relatively 

large scatter among the test data or the fact that the value of tj b given 
by equation (10) is approximate for tubes of moderate length. 

Relatively large scatter can also be observed for the 606I-T6 
aluminum-alloy test data of Moore and Holt (ref. 15; also listed in 
table 3 and shown in fig. 3)- These tests were conducted on both long 
tubes and tubes of moderate length and are so designated in table 3 and 
figure 3* The theoretical values of t^ s were computed from a typical 

stre s s - strain curve for 606I-T6 aluminum alloy given in reference 16 and 
having the same yield properties as listed in reference 15. In this 
case, the test data again favor the use of the maximum-shear law to 
transform the axial data with approximately equal scatter of the test 
points about the theoretical line. 


DISCUSSION 


In discussing the correlation between the theoretical plasticity- 
reduction factors and the available test data, it is convenient to sum- 
marize the results as shown in table 4. 


Compressive Buckling 

For the limited amount of test data on 2017-T4 and 606I-T6 aluminum- 
alloy tubes in compression, it appears that the plasticity-reduction 
factor given by equation (4) is in substantially good agreement with test 
data. Considerably more weight must be placed on the 2017-T4 data as 
compared with the 606I-T 6 data since compression stress-strain data were 
given for the former, whereas such data for the latter were estimated 
from other sources. Thus, the agreement obtained for the 2017-T4 data 
can be interpreted as excellent support for the theoretical value of T] c . 

As indicated previously, the experimental values of t] c for the 
test data were computed using equation (3) which is based on classical 
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stability concepts. Since test data, on elastic buckling of compressed 
cylinders generally fall considerably below this theoretical value, it 
appears fruitful to discuss the implications of the apparent agreement 
of test data on cylinders which buckle inelastically with equation ( 3 ) 
after inelastic buckling effects have been accounted for by use of tj c . 

The cylinders used for the inelastic tests were drawn seamless tubes 
with R/t values less than 50. Consequently, these tubes probably .con- 
tained very small geometrical imperfections and were relatively free of 
residual stresses. By contrast, tests on elastic buckling were generally 
conducted on cylinders fabricated from flat sheets with R/t values 
ranging from 200 to 3*000. Many of these tests were conducted on cyl- 
inders of very thin sheet stock and therefore the geometrical imperfec- 
tions would be expected to be very much greater than those in drawn tubes. 
Thus, it is probably safe to conclude that the initial imperfections for 
the inelastic cylinders were considerably less than those for the elastic 
c ylin ders. Provided the plasticity-reduction factor is correct, the ini- 
tial imperfections were apparently of such a small magnitude for the 
inelastic cylinders that the buckling stress is adequately predicted by 
classical small-deflection theory. 

Thi s apparent agreement with classical theory for inelastic cylin- 
ders may have some bearing on the currently held views concerning the 
lack of agreement of test data on elastic cylinders with Hmall-def lection 
theory. 

According to the energy criterion of buckling used by TBien (ref. 21) 
for perfect elastic cylinders in a rigid screw-powered testing machine 
(the condition leading to the highest buckling stress) 


a cr = 0.37Et/R (25) 

It is the contention of this theory that the small amount of energy nec- 
essary to tri gg er the jump to large deflections is available in the vibra- 
tions of the testing machine, for example. 

Donnell, a n d Wan (ref. 22) have maintained that the presence of geo- 
metrical imperfections and residual stresses rounds off the sharp peak 
in the stress end- shortening curve of large -deflection theory and there- 
fore failure is observed at loads considerably below the classical value. 

It is essential to realize that the highest buckling stress of the 
energy criterion is given by equation ( 25 ) and therefore this theory does 
not admit the possible realization of a buckling stress as high as that 
given by equation ( 3 ). On the other hand, as the imperfections become 
very small, equation (3) is approached as a limit in Donnell’s interpre- 
tation. Since equation (3) was apparently confirmed by the test data 
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shown in figure 1, for which the imperfections were probably very small, 
it would appear that Donnell's interpretation is definitely favored over 
the energy criterion in this case. 


Torsional Buckling 

Of the test data available to evaluate the theoretical plasticity- 
reduction factor for torsional buckling given in figures 2 and 3 and 
summarized in table h } it would appear that considerable weight should 
be placed on the 2017-T4 data. This .is due to the fact that a wide range 
of test data as well as shear stress-strain data was given. For the 
chrome -molybdenum tubes only a small amount of test data could be used 
because the yield stresses and, stress-strain properties varied consider- 
ably from those of the shear stress-strain curve given. Therefore, it 
would, appear that excellent confirmation of t) s was obtained for the 

2017-T4 data, whereas the chrome -molybdenum data were too few and too 
variable in stress-strain characteristics to permit any definite con- 
clusions to be drawn for this material. 

The 605I-T6 and 606I-T6 data of figure 3 are useful in providing a 
means of checking the theoretical value of T) s utilizing the maximum- 

shear and octahedral-shear laws to transform the axial stress-strain 
data to shear stress-strain data. Although the conclusions to be drawn 
are handicapped by relatively wide scatter of test data and by the nature 
of the stress-strain data available for correlation (table 4), it would 
appear that the data are in better agreement with the use of the maximum- 
shear law than with that of the octahedral-Bhear law to transform the 
axial data. 


-SUMMARY OF KESUI1S 


The following conclusions were derived from, a theoretical and experi- 
mental investigation of the compressive and torsional buckling of thin- 
wall cylinders in the yield region: 

1. A general set of equilibrium differential equations for plastic 
buckling of circular cylinders has been derived based on deformation 
stress-strain relations and the no-strain-reversal buckling model. Fur- 
thermore, an attempt has been made to remove a difficulty associated with 
using the no-strain-reversal model in conjunction with classical stability 
concepts . 

2. The plasticity-reduction factors for inelastic buckling of long 
cylinders under compressive or torsional loadings have been derived. 
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It is shown that these factors are in satisfactory agreement with test 
data when satisfactory compression and shear stress-strain data are 
available . 

3- Both the maximum- shear an d octahedral-shear plasticity laws were 
vised in transforming axial stress-strain data to shear stress-strain 
data for torsional buckling of cylinders in conjunction with a theoreti- 
cally derived plasticity-reduction factor based on the octahedral-shear 
law. In such cases, results obtained by use of the transformed shear 
data based on the maximum-shear law are in better agreement with test 
data than those based on the use of the octahedral-shear law. 

4. Compression test data on tubes which probably contained small 
geometrical imperfections correlated very well with the critical stress 
predicted by classical small-deflection theory after the theoretical 
correction for inelastic buckling had been incorporated. This correla- 
tion is viewed as a factor favoring the imperfection interpretation 
(Donnell ) of tests on elastic cylinders over the energy interpretation 
(Tsien) . 


Research Division, College of Engineering, 
New York University, 

New York, N. Y., October 14, 195^- 
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APPENDIX A 

PLASTICITY CONSIDERATIONS 


In the following derivations, assumptions have been employed which 
appear to have resulted in the best agreement between theory and test 
data on inelastic buckling of flat plates with various geometrical bound- 
ary conditions and, types of loading. 


A fundamental hypothesis of plasticity -theory is that the stress 
intensity (Xj_ is a uniquely defined, single-valued function of the strain 

intensity ei for a given material when the stress intensity increases 
(loading) and is elastic when it decreases (unloading). The definitions 
of the. stress and strain intensities theoretically can be chosen from a 
manifold of rotational 1y invariant functions. Two such functions, the 
maximum-shear and octahedral-shear laws, have been useful. 


For the octahedral-shear law, the stress and strain intensities 
can be defined as follows: 



(Al) 


e i 



+ e. 


e x e y + 



(A2) 


With the assumption that the principal axes of stress and strain 
coincide, the secant modulus can be defined as 

E s = °i/ e i (A3) 

Furthermore, by use of deformation-type stress-strain laws together with 
the assumption of plastic isotropy and the idealization that Poisson’s 
ratio is equal to l/2 for both the elastic and the plastic region, the 
following simplified two-dimensional stress-strain laws are obtained: 


= 


“ Cl/*)*] 


(A4) 
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E^[°y “ (1/2)*] 

(A5) 

n 

"S' 

tn 

(A6) 


Inelastic-Buckling Considerations 

All of the foregoing assumptions form the basis for solution of 
plasticity problems in general. For the specific problem of ine lastic 
buckling, it is necessary to make an additional assumption concerning 
the stress distribution at the instant of buckling. 

From the standpoint of classical stability theory, the equilibrium 
differential equations are formulated on the basis that at the buckling 
load an exchange of stable equilibrium configuration occurs between the 
straight form and the slightly bent form. Since the load remains con- 
stant during this exchange, a strain reversal must occur on the convex 
side, and, therefore, the buckling model leading to the reduced-modulus 
concept for columns is correct theoretically. 

Practical columns and plates invariably, contain initial imperfec- 
tions and therefore axial loading and bending proceed simultaneously. 
Since in the presence of relatively large axial compressive stresses 
the bending stresses are generally small, no strain reversal would be 
expected to occur and the incremental bending stresses in the inelastic 
range are given by the tangent-modulus model. However, the bent form 
is the only stable configuration in this case and therefore use of equi- 
librium equations based on perfect columns, plates, or shells is clearly 
unjustified. 

Partially to remove this difficulty, Stowell has assumed that the 
straight form of the plate or column is stable until buckling occurs 
(ref. 5). At buckling, infinitesimal bending is assumed to proceed 
simultaneously with a corresponding infinitesimal increase in b.y1b 1 
loading so that the plate is not subjected to a strain reversal and 
remains inelastic. Again this model poses an essential difficulty since 
classical stability theory is based on the assumption that the a-xiai 
loading remains constant during the buckling process. 

In appendix B, in which the equilibrium equations are considered, 
an attempt is made to remove this difficulty by showing that the infin- 
itesimal increase in load associated with the no -strain-reversal model 
contributes higher order terms than those generally considered in the 
equilibrium equation. This is by virtue of the fact that the axial loads 
are multiplied by first or second derivatives of the displacements and 
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therefore products of the Incremental load increase and these deriva- 
tives result in second-order terms. 


Incremental Forces and M o ments 

When buckling occurs, the displacements vary slightly from their 
values before bu ckling . The resulting strain variations arise partly 
from variations of middle -surface strains and partly because of bending 
strains. These result ing variations of stresses have been considered 
by Ilyushin (ref. 3) and StoweH (ref. 5). Using the assumption that 
no part of the plate is unloaded, Stowell has derived the variations 
of the mome nts durin g the bu ckling process. The variations of the middle - 
surface forces can be derived directly from this work. 

When the variations of the forces and moments are denoted by 
primes ( 1 ), the followin g relations apply to fully plastic plates during 
buckling: 


V - b[Ai€i + ( 1 / 2^62 - (l/2)A 15 e^ 

(AT) 

My* = bJ^ 2 6 2 + (l/2)A/y| e-| - (l^jA^e^J 

(A8) 

N xy' * |^3 e 3 ” C 1 / 2 )^ 6 !. " ( 1 /2)A 52 e^] 

(A9) 

Mjj.' = + (l/2)Ai pXp - (l/2)A]_^^| 

(A10-) 

My 1 = -d|a 2 X 2 + (l/2)A2iXi - (l/2)A 25 X^J 

(All) 

Mxy' = -§^ 5 X 3 - (l/2)A 51 Xi - (1/2)A 52 X2] 

(A12) 

In equations (AT) to (A12), e-j_ and e 2 8X6 middle-surface normal 
strain variations and is the middle-surface shear strain variation; 

and Xg> are the changes in curvature and X 3 is the change in 
twist. Furthermore, the plasticity coefficients are defined as follows: 
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A 1 = 1 ** («*s 2 A) 

*2 - 1 “ ( aa y 2 / ]| ') 

Aj = 1 - err 2 

a 21 = a 12 = 1 “ (o^x 0 ^/ 2 ) 


(A13) 


where 


a 31 ■ a 13 ■ aa x T 
a 32 = ^3 = “V 

“ = (sMf 1 - (*»/*)] 


The axial rigidity is 

B = kB s t/3 (Alii-) 

The bending rigidity is 

D = E s ty 9 (A15) 

In the elastic region, a = 0 and, therefore, A i = A 2 = Aj = 
a 12 = 1 a-nd a 13 = a 23 = O* By replacing the definitions of equa- 
tions (AllO and (A15) which are for a fully plastic plate by 
B = Ety^l - v e 2 ) and D = Wc^jl2.{l. - v e 2 ), respectively, and replacing 
the coefficient (l/2) by v e , equations (A7) to (A12) reduce to the 
familiar relations for the elastic plate. 
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APPENDIX B 

EQUILIBRIUM C0NSH2ERATI0NS 
Elastic Buckling 


Donnell's equations (refs. 13 and 18) for elastic "buckling of t hin - 
wall circular cylinders have been used with, a considerable degree of 
success in buckling problems . Therefore, in this investigation of ine- 
lastic buckling of long circular cylinders under compressive and torsional 
loads, an extension of Donnell's equations is considered. 

The middle-surface strain variations and curvature changes that 
occur during buckling of a circular cylinder are related to the dis- 
placements as follows: 

e-j_ = 3u/3x 1 


e 2 = (3 v/R 30) + (w/R) 



X x ="3^w/3x^ 


„ >2 L 2^ n 2 

*2 = 0 w/R 39 

3x 30 


The following simplified equilibrium equations as derived by Donnell 
(ref. 13 ) neglect certain terms which are of small magnitude when the 
circular cross section of the cylinder is distorted during buckling. In 
cases in which the ctosb section retains its circular shape during buck- 
ling the neglected terms are generally of same importance. 


ZF X 


ay ay 

3x R 30 


0 


EF e 


ay 

R 39 


3N- 


xy 


3x 


(B2) 


+ 


0 


(B3) 
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ZF z 


aV ; aV + H afs +2N a 2 * 

dx 2 R 5x 50 E 2a 0 2 R X Sx 2 ^ R 5x 50 



In equation (b 4), N x , Nxy, Ny, and p are prescribed external 
loadings, with N x and Ny positive in compression. The terms con- 
taining a prime are the variations associated with buckling. 


Inelastic Buckling 

In appendix A, it was indicated that use of the no-strain-re verBal 
model for inelastic buckling poses the difficulty that the external load 
must increase slightly during the buckling process. This is at variance 
with classical stability concepts which require the external load to 
remain constant during buckling. Therefore, some justification is nec- 
essary in order to use equilibrium equations based on constant external 
loads for inelastic buckling problems in which the external load must 
increase slightly. An attempt has been made to clarify this point which 
has been overlooked by previous investigators. 

It is assumed that, in the inelastic-buckling process, the external 
loads increase slightly. Denoting this increment by SN, the external 
loads are increased as follows: N x + 5N X ,. + SN^y, N y + 8My, and 

p + Sp. The terms containing a prime in equations (B2) to (b 4) are the 
middle -surface force and bending -moment variations arising from bending 
and twisting of the plate at buckling. Therefore, the slight increase 
in exte rn al load represented by 8N can have only a negligible influence 
upon the primed terms in equations (B2) to (B4). 

In equation (B^), the external loads R x , N^y, Hy, and p appear. 
If these loads are replaced by R x + 5N X , . . ., p + 5p, then terms 

such as 6 N x ^5^w/5x 2 ^ and 5p appear which are clearly of higher order 

than those terms appearing in equation (b 4) and can be neglected. Thus, 
it appears permissible to conclude that the slight increase in load 
required for the plate, to remain inelastic during buckling is compatible 
with the use of equilibrium equations based on classical stability 
concepts . 
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Equilibrium. Equations 

By use of the force and moment variations (eqs. (A7) to (A12)) and 
the strain-displacement relations (eqs. (Bl)) the equilibrium relations 
(eqs. (B2) to (Bit)) can be written in terms of the displacements u, v, 
and w and their derivatives: 


. Z 2 * a 15 3^ . a 3 3^ a 13 3^ f A 12 k 3\ 3^ 

1 Sx 2 2 E 5x S0 4 E 2 3e 2 ^ 3X 2 \ 2 4/R 3x '30 


a 23 3 ^v ^L2 3w _ ^2£ 3w 

4 R^G 2 2 R 3x 4 R 2 30 


(B5) 


3 2 y k 23 3^v 3^v ^S.3 3 2 u 

2 R 2 30 2 2 R 3x 30 4 ^x 2 it ^x 2 

A 23 3^ a 3w A 23 3w _ Q 
It r 2 3 0 2 2 R 2 30 " k R 3x 


\2 3 2 u 

^ 2 4 /R 3x 30 


(B 6) 


3**w 

3x^ 


3^w 


a*Hr 


Al ^ - A » + ^ + A3 Wa 


-a 23 , S V . ■ + 
P r53x 3©3 


A2 


3\ 


R^O 1 *- 


/ A 12 3u 
A 2 3x 


^23 3u ^3 3v 3v 

— I ■ - -v 1 + A2 " AT T- + A2 “ I + 

If Bd9 4 3x R 30 


0 


3^ 


N x^ +2N ^r^ Ti ’y^ 2 


3 2 v f . „ T 3^ , „ 0 

+ Ky - 6 --— + P = 0 


(B7) 


Equations (B5) to (B7) constitute a basic set of equilibrium dif- 
ferential equations for plastic buckling of circular cylinders. In 
the elastic case, A^_ = Aq ~ A 3 = A 12 = an ^- A 13 = ^23 = and, by 
properly accounting for Poisson’s ratio, equations (B5) to (B7) reduce 
to the following: 
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~ v e + v e . v e &tr _ Q 

Sx 2 2 r 2^ 0 2 2 H 9* de H &c 


(B8) 


~ v e a^v + v e a^u aw 

R 2 ae 2 2 ax 2 2 R R 2 ae 


= 0 


(B9) 


D A + lfv e |a + J5- + 2\ + H x & + a &- + H, 

R\ e ax r ae Ry x ^ ^ * ax ae 


a 2 * 


y ^ +s=o 


(BIO) 


By suitable manipulation of equations (b 8) to (Bio), the above set 
of equations can he reduced to a single equation in deflection w known 
as Donnell’s equation: 



2Nxy 


a 2 * 

r ax ae 


+ % 


a 2 ^ 

R 2 ae 2 



■ 0 


(bh) 
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APPENDIX C 

AXIAL COMPRESSIVE BUCKLING OF A LONG CYLINDER 


For a long circ ular cylinder subjected to axial compression, 

Oy = t =0. The value of the term which appears in equations (A13) 

is given by equation (Al) as discussed in reference 5 and, therefore, 
for th is case = a . K . Thus, the plasticity coefficients reduce to 



1 ?£ 


k E s 


Ag = A^ = A-| p = 1 
A 13 = A 23 “ 0 


(Cl) 


Consequently, equilibrium equations (B5) to (B 7 ) reduce to the following 
equations: 


d 2 * , 

1 c&i 

3 

. 1 J*L- ■ 

dx 2 

* R 2 ^© 2 

+ 4 R &c Se 

2RSx 



3 c&i 

l ^ - 0 

R 2 d0 2 

*3*2 

4 R dx d0 

R 2 B0 


(C2) 


(C3) 



2 d V \ + B/l du + dv + w\ 

R2dx2ae2 R^0V R\2 Sx R S0 Ry 


+ N X 


d^w 

c>x 2 


= 0 


(c4) 


Compressive tests of cylinders which buckle plastically indicate 
that an axi symmetric form of buckling often occurs . Therefore, if this 
mode of instability is assumed, the displacements are inde pendent of the 
0 coor dinat e and all derivatives containing 0 vanish. Thus, equa- 
tions (C2) to (c4) reduce to the following forms: 
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A, 3 “ + 1 ^ -0 

(C5) 

1 Sx 2 2 R 8x 

a/ 

\jo 

II 

0 

(C6) 

dx 2 

Mj. ^ + §4 + l) + N x 2-2 = 0 

3? K V 2 S E / dx 2 

(C7) 


By perfo rming the operation 8/8x on equation (C7) and then using equa- 
tion (C5), a single equilibrium equation in w is obtained: 


da-, ^ + -B-/i - + n x ^ = 0 


(C8) 


A solution to equation (C8) can be written in the following form: 


w = w m sin(i tx/A) 


(C9) 


where 


A = l/m 


Upon substituting the appropriate derivatives of equation (C9) into 
equation (C8) and using the definitions of D, B, and given by 

equations (A15), (Aik-), and (Cl), respectively, and the relation 
N x = c cr t, the following nontrivial solution is obtained: 


cr 


9 \h 4E s ^2 


E t 


2 


9 


(CIO) 
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The cylinder is considered to he long, so that many waves form 
along the length and therefore can be considered as a continuous 

function of Tv. By minimizing equation (CIO) with respect to the wave 
length, 

- f g(j£) ' t™) 


The corresponding half wave length of the buckles 

-HiiW 

For the elastic case, the corresponding solutions are: 


(C12) 


"or * [3(1 

Tv = jt(Bt) 1 / 2 



(C13) 

(cuA) 


By comparing the coefficients which appear in equations (Cll) 
and (C13), it can be observed that the coefficient in equation (Cll) 
can be obtained by substituting a value of l/2 for v e in equation (C13). 

Thus, the following relations can be written which are exact in the 
elastic and fully plastic ranges and result in an excellent degree of 
approximation in the inelastic range; 



(C15) 


(C16) 
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APPENDIX D 

TORSIONAL BUCKLING OF A LONG CYLINDER 


For a long circular cylinder subjected to torsional moments at the 
ends, ct x = Oy = 0. The value of a-j_ which appears in equations (A13) 

is given hy equation (Al) . For this case. 


a ± = (3) 1>/2 t 

The plasticity coefficients reduce to 

A 1 = A 2 = a 12 = 1 

A 13 = A 23 = 0 


(Dl) 


(D2) 


A 3 = Et/ E s 


(D3) 


Consequently, the equilibrium equations (eqB. (B5) to (BY) ) reduce to 
the following expressions: 


d^u a 3 bfci (2 + A^ ^2^ 1 Svf 

dx2 ” 5 * r2^q 2 4 R dx 59 + 2 R dx 


0*) 


+ ^3 + ( 2 * a 3 ) d^u + dw _ 

p2^02 4 ^2 4 R Sx S0 R d0 


+ (l + A 3 )--^£ + 

\bx k ' '''R^x 2 ^ 2 R^0^ 


I . b/i Su , dv , w\ 
n R\2 3x RSfl R/ 


__ 

SNyy v v — 0 

** R Sx d0 


(D5) 


(D6) 
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By multiplying equation (d 6) by R 2 /b and letting 


D _ t £ 


BR 2 12R 2 


_ N xy _ T cr 
9 B (k/3)B e 


(D7) 


Equation (l>6) can be reduced to the more convenient form: 


P 


■o 1 !- d**w . ft . ft \t> 2 d 4 w . d 4 w 
R + (1 + A 3 jR 


+ S fa + |z + W + 2cpR J&- = 0 
2 Sx de ^ Sk ae 


(B8) 


Following the method of solution for torsional elastic buckling of 
a cylinder as given by Timoshenko - (ref . 1), the following relations are 
used for the displacements: 


u = %n cos! 


<T ' ” 6 ) 

<T - * 9 ) 

w = Wmn sin0£ - noj 


v = Yum cos 


(D9) 


The angle which the helical, buckle makes with the original generator of 
the cylinder is given by 


tan a) = ~h/n (DIO) 

Since the cylinder is assumed to be long, the boundary conditions at the 
ends are relatively unimportant and, therefore, equations (D9) can be 
used although they do not satisfy the usual boundary conditions of simple 
support or clamping. 
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The substitution of the appropriate derivatives of equations (D9) 
into equations (D4-), (D5), and. (d 8) results in the following matrix 


( 

«5fl 

-1 

CM 

£ 

1 

( 2 + A 3) m 



V T k 


2 A 


( 2 + ^An 

in 2 + tl A 

— n 

V™. 

V if / 

V n + k h J 


v mn 


-n 

- [l + 2cpAn + PA* 4- + 
2p^l + A^A^ 2 + pn*j 



A nontrivial solution requires the determinant to vanish, independently. 
Upon expanding the determinant, the following result is obtained: 

A* 1- + p|n^ + ^1 + Ajj'hPa 2 + A**| + A^A^j 

cp = j= = * (Dll) 

2AnjA^n^ + (3 - A^A^ 2 + AjaM 

As in the elastic case, very small values of A yield the smallest 
value of cp. Furthermore, by assuming p to be small, equation (Dll) 
simplifies to 

„ _ (?A)^ * fra pi 


Note that the term Aj cancels out of equation (D12) u n d e r these 

assumptions and, therefore, the solution for cp is independent of the 
value of Aj given by equation (D3). 

For a long cylinder a minimum value of cp is obta i n e d for n = 2 
which corresponds to the familiar two-lobe buckling. Therefore, equa- 
tion (D12) becomes 


(3A)A 4 + 256p 


(D13) 
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By differentiating with respect to X, 



(raA) 


Upon substituting equation (D]A) into equation (D13) and simplifying. 


* - 


(D15) 


By use of equations (D 7 ) 


T cr = |(3)“ 1 ^ s (t/R ) 5/2 


= 0.338E s (t/R) 


5/2 


(Dl 6 ) 


The corresponding elastic solution can be obtained by substituting 
1 - v e 2 for the coefficient 3A in. equations (D13) and (D 7 ): 



(D17) 


(D18) 
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By Identifying the coefficient hfe which appears in equations (D13) 

and (D7) as the limiting value of (l - v^) with v = l/2, it is 
then possible to combine equation (Dl6) with (D17) and equation (Dl4) 
with (Dl8): 


= 0 . 272(1 - v 2 ) 3// \ s (t/R) 3 / 2 


L cr 


"K = 



(D19) 


(D20) 


Equations (D19) and (D20) are exact for the elastic and plastic ranges 
and constitute excellent approximations in the inelastic range. 
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TABLE 2 

TEST DATA ON TORSIONAL STRENGTH OF LONG TUBES FROM REFERENCE 19 
(a) 2017-T4 aluminum-alloy tubes 


1 / a 

t/a 

Tf, ksi 

% 

19.9 

0.02235 

23.1 

0.792 

19-9 

.02516 

23.4 

.69 8 

60.0 

.02848 

24.2 

.580 

19.9 

.03231 

24.9 

.491 

59.8 

.03242 

26.0 

.513 

39.9 

.02195 

22.8 

.830 

13.3 

.02903 

25.1 

•596 

39-9 

.02896 

24.1 

.572 

13.3 

.03273 

26.4 

.512 

4o.o 

.03312 

26.1 

.506 

30.0 

.03350 

25.7 

.491 


Cb) Chrome -molybdenum, tubes 


z/a 

t/a 

T yieia> kBi 

ksi 

^s 

25.3 

0.04055 

48.6 

50.6 

0.262 

25.3 

.o4o4o 

47.9 

50.4 

.261 

79-9 

.04020 

49.5 

51.1 

.264 

30.1 

.02330 

47.2 

43.8 

.530 

12.6 

.03486 

49.8 

47.4 

.309 
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TABLE 3 


TEST BATA ON TORSIONAL STRENGTH OF TUBES 
(a) 6051 -T6 aluminum-alloy tubes of moderate length from reference 20 



d/t = 

77 

d/t = 

10lf 

d/t = 

139 

1 / d 

r - f , ksi 


Xf , ksi 

ns 

Xf , ksi 

n B 

1 

21.6 

0.412 

18.9 

0.525 

18. If 

0.728 

1 

21.9 

.If l8 

19.5 

• 5!f2 

18. If 

.728 

2 

21. If 

.580 

18.8 

.7l|0 

17.0 

.950 

2 

21.9 

.594 

19.7 

.765 

16.5 

.923 

4 

21.2 

.810 

17.7 

.98lf 

12.6 

1.00 

4 

21.3 

.8llf 

18.6 

1.03 

12. If 

.982 

8 

20.1 

1.09 

13.9 

I.09 

9.2 

1.035 

8 

19-2 

l.Olf 

llf.O 

1.10 

9.2 

1.035 

16 

13.5 

1.03 

9-8 

1.09 

6.5 

1.02 

16 

12.5 

•957 

9.6 

1.08 

6.5 

1.02 


(b) 606I-T6 aluminum-alloy tubes from reference 15 


Length 

(a) 

d/t 

Z/d 

Xf, ksi 

n B 

M 

80.6 

5.7 

9 

0.816 

M 

80.6 

17.0 

W&Km 

.813 

L 

80.6 

27.6 


•795 

M 

80.6 

11.5 

• 

.987 

M 

58.8 

5-7 

18.6 

. 5^3 

L 

58.8 

17.0 

17. 4 

.927 

L 

58.8 

27.6 

15.2 

.810 

M 

39-4 

5.7 

22.7 

.403 

L 

39 . If 

17.0 

22.3 

.648 

L 

39 A 

27.6 

22.0 

.640 

M 

60.6 

H .5 

20.2 

.928 


a M, moderate length; L, long. 
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Material 

TnnSIng 

Figure 


T1 

Stress-strain 

data 

Eemarks 

2017-T4 

Compression 

1 

Eq. 

(M 

Compression 

Excellent agreement; 
limited test data 

666L-T6 

Compression 

1 

Eq. 

00. 

a Urt-ir 

cy 

Fair agreement; limited 
test data 

2017-T4 

Torsion 

2 

Eq. 

(17.) 

Shear 

Excellent agreement; 
vide range of test data 

Chrome - 
molybdenum 

Torsion 

2 

Eq. 

(17) 

Shear 

Theory high; limited test 
data 

6051-T6 

Torsion 

3 

Eq. 

(10) 

Tension 

Fair agreement; moderate 
length tubes 

6061-T6 

Torsion 

3 

Eq. 

(10) 

a °y 

Good agreement; moderate 
ftnrl long tubes 


B Tensile and compressive yield stress given, 
ref. 1 6 for these yield values. 


Stress-strain data obtained from. 
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Figure 2.- Plasticity-reduction factoi 
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Figure 3>- Plasticity-reduction factors for torsional bu ckling of tubes. 
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